An atom interferometer enabled by spontaneous decay 
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We investigate the question whether Michelson type interferometry is possible if the role of the 
beam splitter is played by a spontaneous process. This question arises from an inspection of trajec- 
tories of atoms bouncing inelastically from an evanescent-wave (E W) mirror. Each final velocity can 
be reached via two possible paths, with a spontaneous Raman transition occurring either during the 
ingoing or the outgoing part of the trajectory. At first sight, one might expect that the spontaneous 
character of the Raman transfer would destroy the coherence and thus the interference. We inves- 
tigated this problem by numerically solving the Schrodinger equation and applying a Monte-Carlo 
wave-function approach. We find interference fringes in velocity space, even when random photon 
recoils are taken into account. 

PACS numbers: 03.65.Yz, 03.75.Dg, 39.20.+q 



I. INTRODUCTION 

Spontaneous emission is generally considered a detri- 
mental effect in atom interferometers. The associated 
random recoil reduces or even completely destroys the 
visibility of the interference fringes. In this paper we 
describe an atom interferometer where the beam splitter 
works by means of a spontaneous Raman transition. Our 
central question will be whether one can observe inter- 
ference in such an interferometer which is enabled by a 
spontaneous process and where decoherence is built into 
the beam splitting process from the start. 

The role of spontaneous emission in (atom) interferom- 
eters has long been connected to the concept of which- 
way information, ultimately tracing back to Heisen- 
berg's uncertainty principle Q. Feynman [j| discussed 
a gedankenexperiment using a Heisenberg microscope to 
determine which slit was taken by a particle in a Young's 
two-slit interferometer. The scattered photon, needed 
to determine the position, spoils the interference pattern 
due to its associated recoil. Similarly, spontaneous emis- 
sion in an (atom) interferometer may provide position 
information on the atom, while at the same time random- 
izing the momentum and spoiling interference 0, Q . On 
the other hand, resonance fluorescence, including that 
from spontaneous Raman transitions, is as coherent as 
the incident light in the limit of low saturation HQ. Ex- 
perimental demonstrations have been given by Eichmann 
et al. |3 and Cline et al. |8j. An experiment by Diirr 
et al. 9] has made it clear that availability of which-way 
information should be conceptually separated from the 
presence of random recoils. Which-way information can 
be obtained without random recoils, nevertheless lead- 
ing to loss of interference. Here we show that random 
recoils do not lead to loss of interference as long as the 
spontaneously emitted photons do not yield which-way 
information. 

Our proposed interferometer is based on cold atoms 
that reflect from an evescent-wave (EW) mirror. Our 
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FIG. 1: Atoms approaching an evanescent- wave mirror in 
state |1) can undergo a spontaneous transition to state |2), ei- 
ther on the ingoing or on the outgoing part of the trajectory. 
These two paths possibly interfere, with a phase difference 
depending on the transition point z. 



model atom is a typical alkali atom with two hyperfine 
ground states. During the reflection from the EW mirror 
the atoms can make a spontaneous Raman transition to 
the other hyperfine state. When the repulsive potential 
experienced by the final state is lower, the atoms lose 
kinetic energy and hence bounce inelastically from the 
potential 0, 0]- This Sisyphus process has been in- 
vestigated previously 01 an d the resulting final velocity 
distribution has been shown to be a caustic [l^. remi- 
niscent of the rainbow. Furthermore it is used in several 
experiments as the loading process of (low-dimensional) 
traps for atoms 0, 0, . Cognet et al. 0] observed 
the analog of Stiickelberg oscillations in the transverse ve- 
locity distribution of atoms that reflect elastically from 
a corrugated EW potential. However, no stochastic or 
incoherent processes were involved. 

In our case the final velocity of an atom depends on the 
position where it made the Raman transfer . Looking 
at the trajectories in detail we see that each final velocity 
can be reached by two trajectories, as is shown in Fig. 
n An atom can be transferred to the second state on 
the ingoing or the outgoing part of its trajectory. Inter- 
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ference will manifest itself in the velocity distribution of 
the reflected atoms, since the outgoing velocity depends 
on the transition point z. Note that the beam splitter, 
its role being played by a spontaneous Raman transition, 
is highly non-unitary: atoms are only transferred from 
state |1) to state |2) and not vice versa. 

This paper is structured as follows. We first present a 
semi-classical picture, and use it to make qualitative pre- 
dictions about the behavior of the interference effects, 
if present. The question whether interference is possi- 
ble will be answered by solving the Schrodinger equation 
in two different ways. The first approach will employ 
stationary analytical solutions of the time-independent 
Schrodinger equation, but is limited to monochromatic 
wave functions. The second approach will propagate 
a wave packet by numerically integrating the time- 
dependent Schrodinger equation with random quantum 
jumps describing the Raman transitions. The last section 
deals with experimental considerations. 

II. SEMI-CLASSICAL DESCRIPTION 

The analysis will be for a two-level atom, and only 
the motion of the atom in the direction along the EW- 
potential gradient will be considered. The calculations 
throughout this paper will assume a low saturation pa- 
rameter, so that depletion of the initial state can be 
neglected. An atom in state |1) that reflects from an 
EW mirror experiences a potential VLexp(— 2kz), with 
kT 1 the decay length of the EW field. The atom's 
trajectory through phase-space is given by z%(v) — 
(— 1/2k) In \ijn/2V\)(v 2 — v 2 )] , with v\ the velocity with 
which the atom enters the potential, see Fig. [3 Af- 
ter a transition to state |2) in point A or A' the atom 
experiences a potential j3Vi exp(— 2kz), with j3 < 1 the 
factor by which the potential energy is reduced after 
the Raman transition. The atom continues its way 
through phase space on a new trajectory z-i (v), given 
by z 2 {v) = (-l/2n)\n[(m/2/3V 1 )(vf - v 2 )], with v { the 
asymptotic velocity with which the atom leaves the po- 
tential. The final velocity can have any value between 
Wf = \/J3vi, which corresponds to a transfer in the turning 
point, and Vf = v- u which corresponds to a transfer out- 
side the EW potential. The final velocity depends on the 
atom's velocity Vt at the moment of the transfer to state 
1 2). This dependence is given by v 2 — v 2 + [3(vf — v 2 ), 
from which it is clear that two values ±u t lead to the same 
final velocity V{. The two transitions lead to two possible 
trajectories through phase space, as shown in Fig. |2{a). 
The phase difference Aip between the two trajectories is 
given by 

m f Vt 

A( P=jr J {z 1 {v)-Z 2 (v))dv, (1) 

This phase difference is proportional to the area between 
the two curves, indicated in gray in Fig. |2{b). From 
the evaluation of the integral for various parameters we 
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FIG. 2: Phase-space trajectories of atoms being repelled by an 
evanescent potential. Atoms initially in state |1) can be trans- 
ferred to state 1 2) and continue on a different path through 
phase space, (a) Depending on the initial shape of the wave 
packet (e.g. the light and dark grey areas) which-way infor- 
mation can be obtained or not. (b) The accumulated phase 
difference between these two paths, indicated by the enclosed 
gray area, may give rise to interference effects, (c) A momen- 
tum kick due to the spontaneous recoil gives rise to an extra 
phase contribution. 



learn that the fringe period decreases for increasing initial 
velocities v- u for increasing final velocities Vf, for smaller 
j3 and for larger decay length (smaller k). 

In a semi-classical picture the atom can be treated as a 
wave packet which is subject to Heisenberg's uncertainty 
relation AzAv z > h/2m. The distribution of uncertainty 
between position z and momentum mv z is determined 
by the experimental preparation procedure of the atoms. 
For a wave packet that initially has a large spread in mo- 
mentum, it is not possible to unambiguously determine 
the phase difference between the two possible paths, since 
the wave packet is spread out over several classical tra- 
jectories through phase space. This is indicated by the 
light grey area in Fig. E[a). It is, however, possible to de- 
termine whether the transfer to state |2) is on the ingoing 
or outgoing path, by observing the timing of the spon- 
taneously emitted photon. Therefore, it is not expected 
that wave packets with this shape show interference. On 
the other hand, a minimum uncertainty wave packet with 
a narrow initial momentum spread will more closely fol- 
low a classical trajectory through phase space. This is 
indicated by the dark grey area in Fig. |2{a). The phase 
difference between the two paths is well defined. The two 
points in phase space, A and A', where a transition to the 
final trajectory is possible are covered simultaneously by 
the wave packet. Thus no which-way information can be 
obtained by observing the time of emission of the spon- 
taneously emitted photon. The initial trade-off between 
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position and momentum uncertainty in a bandwidth lim- 
ited wave packet determines whether interference can be 
a priori excluded or not. 

The random direction of the spontaneously emitted 
photon can be taken into account in the motion of the 
atom by a random momentum jump. This makes the 
atom propagate on a different trajectory through phase 
space than it would have without the random recoil. The 
momentum changes are indicated by horizontal arrows in 
the phase-space diagram of Fig. Hfc). For a single atom, 
or a collection of distinguishable atoms, the spontaneous 
recoil could be measured by detecting the direction in 
which the photon was emitted. Due to this possibility 
there will be a set of interference patterns, one for each 
recoil direction. By disregarding the information present 
in the scattered photons, we probe the incoherent sum of 
all these interference patterns. The phase difference be- 
tween the two paths is different with respect to the recoil 
free case, and it depends on the direction of the recoil. 
It is indicated by the gray areas in Fig. Hfc). In order 
for the interference to be experimentally observable the 
difference between the interference patterns with a cer- 
tain recoil direction should not be too large. This means 
that the phase difference between recoil components in 
the ±z directions should be less than it. For larger final 
velocities these phase corrections get larger as is apparent 
from comparing the areas around the point B' with the 
areas around point A' in Fig. |2Ic). We thus expect the 
visibility of the interference to decrease for larger final 
velocities. 

Note that which-way information can neither be re- 
trieved from a measurement of the frequency uj of the 
emitted photon. This frequency is determined by en- 
ergy conservation: hui = hm(vf — vf) + Huiew — hAi2, 
where wew is the frequency of the evanescent photon and 
hAi2 is the energy difference between states |1) and |2). 
Because the initial and final kinetic energies are equal 
for both trajectories, the frequency of the spontaneously 
emitted photon is equal for both interfering paths. 



III. TIME-INDEPENDENT APPROACH 

The question whether interference is visible will be an- 
swered in this section by considering analytical solutions 
of the time-independent Schrodinger equation 

+ v ie -^Mz) = p-M*), (2) 

2m oz z 2m 

describing stationary states with a total energy p§ /2m on 
the potential V\ exp(— 2kz). It thus describes particles 
with momenta ±p in the asymptotic limit of large z. 
This is one of the few examples where the eigenfunctions 
of the Schrodinger equation are known analytically. The 



solutions are given by 




where K a (-) is the Bessel-K function of order a jl8| . 
These functions are normalized such that the asymptotic 
density is independent of po- They are also given by 0] 
where a different normalization is used. 



This wave function describes the atoms that are inci- 
dent in state |1). After the spontaneous Raman tran- 
sition to state 1 2), the atoms are described by one of 
the eigenfunctions il>2,p(z) with final momentum p in the 
potential V% exp(— 2nz). The final wave function in mo- 
mentum space is given by the overlap integral 

/oo 
^)e- KZ - ikz r2, P (z)dz, (4) 
-00 

where the recoil due to the absorbed evanescent photon is 
taken into account by the factor exp(— kz) and the recoil 
due to the spontaneously emitted photon by the factor 
exp(— ikz), with hk the momentum component of the re- 
coil in the z direction. As already discussed there will be 
interference patterns |</>fc(p)| 2 for every value hk of the re- 
coil. A measurement that disregards the emitted photon 
yields the sum of all these possible interference patterns. 
In this derivation we will assume an isotropic distribution 
of the recoil momentum hk. This is an approximation, 
since the distribution depends on the polarization of the 
spontaneously emitted photon. We will come back to this 
point in section [V] This leads to 

\HP)\ 2 = f +k ° \Mp)\ 2 dk, (5) 
J-k 

with hko the total recoil momentum. 

Fig. |3 shows the behavior of the momentum distribu- 
tion for various values of the z component of the photon 
recoil hk. It is calculated using Eq. |0J, with an initial 
momentum po = 2hko, a potential steepness k = /co/8, 
and potential reduction [3 = 0.2. The main features of 
this figure can be understood from semiclassical argu- 
ments. 

The dashed curved and straight lines demarcate three 
separate classically allowed regions, taking into account 
the recoil. The dashed curve on the left indicates the 
lowest classically reachable momentum, given by p — 
- (hk/p Q ) 2 /(l - To the right of this curve, 
we clearly see a region with interference. In addition we 
see two triangular regions with no interference, demar- 
cated by two straight lines described by p = po ± hk. The 
triangular regions can be reached by a single phase space 
trajectory only. The upper (lower) triangle corresponds 
to a trajectory where the Raman transfer took place on 
the outgoing (incoming) branch. Only the region on the 
left of the triangles is reachable by two trajectories and 
thus shows interference. 

The left dashed curve is reached by atoms that scat- 
ter a photon near the turning point. Note that in this 
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FIG. 3: The behavior of the interference pattern |</>fe(p)| 2 ver- 
sus the final momentum p, expressed in units of the photon 
recoil hko, for various values of the z component of the photon 
recoil hk for parameters po = 2hko, k = feo/8, and /3 = 0.2. 
The dashed curve indicates the classically lowest reachable 
momentum. Interference is visible in the area that can be 
reached by two trajectories, enclosed by the dashed curve and 
straight lines. The triangular regions can only be reached by 
one trajectory and hence no interference is visible. 

case the photon recoil has only a small influence on the 
final momentum p. The final momentum is mainly deter- 
mined by the potential energy near the turning point that 
is converted to kinetic energy. The amount of kinetic en- 
ergy that can be added or removed by the photon recoil 
near the turning point is small because the atomic veloc- 
ity is small. As a result we see only a slight curvature as 
a function of the recoil k. For larger values of the initial 
momentum pq or the ratio (3 = V2 /V\ the left curve will 
become more and more straight. 

As expected, the main part of the momentum distri- 
bution is in the classically allowed regions. The distri- 
butions peak near the lower classical limit (the dashed 
curve), resembling the caustic distribution |13| . For ev- 
ery recoil direction interference is visible. Although the 
region with interference is smaller for larger values of the 
recoil, the remaining interference fringes are present at 
more or less the same final momenta. This indicates that 
the spontaneous recoil does not completely wash out the 
interference. The behavior of the interference does not 
depend on the sign of the recoil, because a photon that 
is emitted on the ingoing part of the trajectory has the 
same effect on the momentum distribution as a photon 
that is emitted in the opposite direction on the outgoing 
part of the trajectory. 

Fig. 01 shows the final momentum distribution, calcu- 
lated by Eq. J5J) for various experimental parameters. 
The results are compared with an evaluation without a 
stochastic contribution. Indeed, the averaging over the 
spontaneous recoil does not destroy the interference pat- 
tern. The small part of the distribution that extends into 
the classically forbidden region is an evanescent matter 
wave. Several of our predictions that were made for the 
general case of a coherent interferometer are noticeable 
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FIG. 4: Final momentum distributions j0(p)| 2 calculated 
for different parameters using the time-independent method. 
Solid lines including the effect of a spontaneous recoil, and 
dashed lines without the effect of recoils, versus the final mo- 
mentum p in units of the photon recoil hko ■ The dotted lines 
indicate the classically allowed region without recoil. Between 
(a) and (b) the initial momentum po is changed, (a) and (c) 
differ in the reduction factor /3, (c)-(f) are a sequence for de- 
creasing decay length k~ . 

in these graphs. Indeed, the fringe spacing decreases for 
larger initial and final momenta, for longer decay lengths 
kT 1 of the evanescent field, and for smaller values of j3. 
Furthermore, as predicted for the case of an incoherent 
interferometer, the visibility of the graphs in which the 
effect of the recoil has been taken into account decreases 
for larger final momenta. 

IV. TIME-DEPENDENT APPROACH 

In the previous section we have shown that the incoher- 
ent nature of the spontaneous Raman transfer does not 
prevent us from observing interference. In this section we 
show that the interference phenomena will also be visible 
for a wave packet with a finite momentum spread. In the 
analysis we closely follow the Monte-Carlo wave-function 
approach poll2ll|. 

We consider the evolution of a diffraction limited wave 
packet in state |1) 

W ^ = 0) = /r^W e,feV ^ (6) 

with initial height zq, initial width a Zl and initial momen- 
tum po = hk z . It is normalized such that J \?p(z, 0)| 2 dz = 
1 and J(z — z ) 2 \ip(z : 0)| 2 dz = a\. The evolution of the 
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wave packet when it reflects from the evanescent-wave 
potential with a potential height V\ at z = is calculated 
by numerically solving the time-dependent Schrodinger 
equation 

ifij^M** t) = ~^^2^ z > *) + Vie~ 3 "*M*> *) (7) 

using the Quantum kernel [22| package in Mathematica 
|23. This results in a wave packet i/^i (2;, t) at time t. 
At a time r a spontaneous Raman transition to state 
1 2) occurs, and the evolution abruptly continues on a 
potential that is a factor (3 lower. Immediately after the 
transfer the wave function in state |2) is described by 

i> r , k {z,t = t) = ^ x {z,T)e- KZ e- ikz , (8) 

where Ji denotes a normalization factor. The two ex- 
ponents are equal to the exponents in Eq. J3J. The 
evolution of this wave function can now be continued up 
to a time t eu &, leading to a wave function Vv,fc(^> *end)- 
When f eru j is large enough the entire wave packet effec- 
tively propagates in free space, so that the momentum 
distribution remains constant. The Fourier transform 

4> T ,k{p) — yWv,fc(-Mend)) (9) 

of the wave packet at this time is the wave function in 
momentum space that endured a momentum kick hk at 
its transfer time r. The transfer rate T(r) at a certain 
time r is given by 

/>OC 

r(r)cx / ^1{z,t)Vx{z)^i{z,t)6z. (10) 
Jo 

We again assume an isotropic distribution of the recoil 
momentum hk. Contributions with different r and k have 
to be summed in an incoherent way. For the momentum 
distribution as a function of the wave vector of the spon- 
taneously emitted photon we get 

/■tend 

\MP)\ 2 « / r(r)|0 Tifc (p)| 2 dr. (11) 
Jo 

A subsequent integration over this wave vector yields 

\HP)\ 2 = [° \Mp)\ 2 dk (12) 
J —ko 

for the momentum distribution of a sample of atoms. 

Fig. [5] shows graphs of |0(p)| 2 for some parameters. 
All calculations are performed with t cn d — 70m /Tik^ for 
which the criterium that the entire wave packet has left 
the potential is fulfilled. Fig. 02a) should be compared 
with Fig. Ufa). The parameters for Figs. Ub)-(d) are 
equal to the parameters used for Fig. 0{c) except for 
the initial width a z of the wave packet and thus the 
momentum spread a p = hja z . Also for a wave packet 
with a finite momentum spread the interference effects 
are present. As expected, the interference fringes are 
more apparent for a wave packet with a smaller momen- 
tum spread. 
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FIG. 5: Interference patterns calculated for different parame- 
ters calculated using the time-dependent method. Solid lines: 
|0(p)| 2 , with the effect of a spontaneous recoil, and dashed 
lines: |0o(p)| 2 , without the effect of a spontaneous recoil, 
versus the final momentum p in units of the photon recoil 
hko. The dotted lines indicate the classically allowed region 
without recoil. Between (a) and (c) the reduction factor f3 is 
changed, (b)-(d) are a sequence where the initial momentum 
spread a p is decreased. 



V. EXPERIMENTAL CONSIDERATIONS 

So far we have considered levels |1) and |2) without 
discussing which physical level they correspond to. In 
reality we usually deal with multi-level atoms, that more- 
over include sub-structure. Each of these (sub-)levels has 
a different interaction with the evanescent field. If more 
(sub-)levels contribute to the signal the predicted inter- 
ference can be washed out. For 87 Rb atoms a convenient 
choice would be the \Fm) — |1, 0) ground state for state 
|1) and the \Fm) = |2,±1) ground states for state |2). 
The evanescent field needs to be linearly (ir) polarized 
and blue detuned with respect to the F = 1 — > F' = 2 
transition of either the D\ or D2 line. Due to selection 
rules, only a transition over the \F'm!) — |2,0) excited 
state contributes to the transition from state |1) to state 
|2). This excited state can decay to either of the |2,±1) 
ground states by emitting a polarized photon. Since 
these states interact identically with the evanescent field, 
their interference patterns will overlap. 

The necessary linear polarization for the EW can only 
be easily obtained using a TE (_L z) polarized incident 
beam. For a circularly polarized photon this distribution 
of wave vectors of the spontaneously emitted photon is 
non-isotropic. The circularly polarized photon has its 
quantization axis along the polarization axis of the EW 
field. The recoil distribution in the z-direction due to 
such a photon is given by ^(3 — (k/k ) 2 ). Intensity dis- 
tributions of dipole radiation are given by e.g. [24|. Since 
this distribution has a maximum for k = 0, for which 
the recoil-dependent interference patterns are most pro- 
nounced as is visible in Fig. |3| the visibility of the inter- 
ference signals will be slightly better than presented in 
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this paper. 

VI. DISCUSSION AND CONCLUSIONS 

The calculations presented in this paper have been per- 
formed for low initial velocities U;. This is because both 
calculation procedures turned out to be limited by com- 
putational resources. For the time-independent approach 
the evaluation in Mathematica [2j| of the Bessel-K func- 
tions of high imaginary order becomes very slow. For the 
time-dependent approach the number of sampling points, 
necessary for the numerical evaluation, becomes too large 
due to the highly oscillatory character of the incident 
wave packet. However, we expect even better signals for 
realistic values for the initial velocity V{, so the calcula- 
tion represents a worst case. 

By numerically solving the Schrodinger equation we 
can now reply with an unambiguous yes to the question 
can the beamsplitter in an atom interferometer work on 
the basis of spontaneous emission? The intuitive objec- 
tions to whether this is possible have been refuted. The 
semi-classical arguments have been confirmed by the full 
quantum- mechanical calculations. Which- way informa- 
tion due to the possibility of detecting the time of emis- 
sion of the spontaneously emitted photon is avoided by 
choosing a sufficiently narrow velocity uncertainty A 



wave packet covers both transfer points in phase space si- 
multaneously if its velocity is defined accurately enough. 
The incoherent nature of a spontaneous emission process 
due to the random recoil direction of the atom is visi- 
ble in all the calculated interference curves, but does not 
lead to a complete scrambling of the interference. For 
larger final velocities, for which the transfer points are 
separated more and the acquired random phase is con- 
sequently larger, the visibility of the interference fringes 
indeed decreases. Furthermore the fringe period indeed 
qualitatively shows the behavior that was predicted on 
the basis of the semi-classical calculations. 

Our analysis also shows that the absence or presence of 
which-way information is not the same as the perturbing 
effect of the recoils due to spontaneous transitions. This 
is most clearly seen in Fig. [21 and is in agreement with 
the viewpoint of Durr et al. [9( . 
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